Excitation of Alfven Waves and Pulsar Radio Emission 

Maxim Lyutikov 
Canadian Institute for Theoretical Astrophysics, St George Street, Toronto, Ontario, M5S 3H8, Canada 

ABSTRACT 

We analyze mechanisms of the excitation of Alfven wave in pulsar magnetospheres as a 
possible source of pulsar radio emission generation. We find that Cherenkov excitation of 
Q^^ , obliquely propagating Alfven waves is inefficient, while excitation at the anomalous cyclotron 

0^ ' resonance by the particles from the primary beam and from the tail of the bulk distribution 

function may have a considerable growth rate. The cyclotron instability on Alfven waves occurs 
Q [ in the kinetic regime still not very closed to the star: r > 50Rns- We also discuss various 

^^ ' mechanisms of conversion of Alfven waves into escaping radiation. Unfortunately, no decisive 

I ■ conclusion about the effectiveness of such conversion can be made. 
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^P ' Interpretations of various observational data tend to place the location of radio emission generation at 

a distance r w 10 — IOORns (e.g., Phillips 1992), though there plenty of claims to the contrary (Kijak et 

Qv ' al. 1999, Smirnova et al. 1996, Gwinn et al. 1997). As was pointed out by Kunzl et al. (1998) and Melrose 

Q^ ■ & Gedalin (1999), even with the most conservative estimates of the efficiency of plasma production in the 

polar caps, the plasma frequency at those heights is much larger that the observed frequency. This argues 
C^h' against radio emission mechanisms that generate Langmiur waves with a frequency near the local plasma 

Q i frequency (Asseo et al. 1990, Whetheral 1997). von Hoensbroech et al. (1998) argued that this implies 

a strong underdense production of particles, but the theoretical foundations of such an assumption are 
^ ' weak. Alternatively, Melrose & Gedalin (1999) argued that in order to restrict emission to small altitudes 

the emission should be generated at frequencies much smaller than the plasma frequency, preferably on 
Alfven waves. They considered excitation of oblique Alfven waves by the Cherenkov resonance with plasma 

p\. • and found that it mostly produces waves with cu ^ ujp and thus cannot resolve the problem. Though 

we agree with their conclusion that Cherenkov excitation of Alfven waves is insignificant in the pulsar 
magnetosphere, we disagree on the reasons why. First, we do not agree with the conclusion of Melrose & 
Gedalin (1999) that Cherenkov resonance of the beam particles with the Alfven waves occurs outside the 
light cylinder for the conventional beam energies. It actually occurs for radii r < 50Rns (see section ||). 
Secondly, they assumed that Cherenkov excitation of Alfven waves occurs in a kinetic regime, while it was 
shown in Lyutikov (1999) that Cherenkov excitation of Alfven waves occurs in a hydrodynamic regime. 

On the other hand, some authors postulated excitation of the Alfven mode and, using the fact 
that Alfven is guided along the curved magnetic field, were able to explain some observational data like 
dependence of the mean profile on the frequency (e.g. Barnard & Arons 1986, McKinnon 1994, Gallant 
1998, Gwinn et al. 1999). The above arguments stimulated us to reconsider the possibility of excitation of 
Alfven waves at the low altitudes in the magnetosphere, including the excitation of waves at anomalous 
cyclotron resonance. 

A major theoretical problem of the theories that produce radio emission on Alfven waves is that Alfven 
waves cannot escape from the plasma and have to be converted into escaping radiation (preferably X mode). 
In section ra we review various possibilities of Alfven waves conversion into escaping modes. 



2. Waves and resonances 

The open field lines of the pulsar magnetosphere are populated by the dense one-dimensional flow 
of electron-positron pair plasma penetrated by highly energetic primary beam with the density equal 
to the Goldreich- Julian density ncj = r2.B/(2 7rec), Lorentz factor 7^ sa 10^ (Arons 1983, Daugherty 
& Harding 1996). The density of the pair plasma is n « ^mtigj ~ lO'^ — lO^ncj, where Am is the 
multiplicity factor which gives the number of pairs produced by each primary particle; its Lorentz factor is 
7p ~ 7fc/AAf = 10 — 10"^. The distribution function of the bulk plasma also has a high energy tail up to to 
the Lorentz factor 74 w 10^. 

Though the pulsar plasma is thought to be relativistically hot, we restrict ourselves to cold plasma case 
which simplifies consideration considerably. Except for the Landau damping of the waves with slow phase 
velocity, thermal effects has only marginal importance for the wave-particles interaction and growth rates 
of instabilities (Lyutikov 1999a). In what follows we consider wave excitation in the plasma frame and then 
use the transformation rules of the Lorentz factors and frequencies from the plasma frame to pulsar frame: 
7' = 27p7, Wp = ^/TpWp, uj'g — ojb, where primes denote quantities in the pulsar frame and 7^ is the Lorentz 
factor of the relative motion of the plasma frame to the pulsar frame. 

In the strongly magnetized pair plasma, in the small frequency limit cj <C tjp, there are two modes 
- transverse extraordinary (X) mode, with the electric vector perpendicular to the k-B plane, and 
quasitransverse Alfven wave with the electric vector in the k-B plane. In the plasma frame the dispersion 
relations of the extraordinary and Alfven modes are (e.g. Lyutikov et al. 1999a). 

fee cos 1 ^ 2 — ' for tj < Wp (1) 
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The plasma normal modes can be excited due to interaction with particles at the Cherenkov 

uj — k\\v\\ = (3) 

and anomalous cyclotron resonances 

UJ -k\\V\\+UB/lres^Q (4) 

where w is a frequency of the wave, k\\ and v\\ are projections of the wave vector and velocity along the 
direction of the magnetic field, ujb is a (nonrelativistic) cyclotron frequency and ^res is a Lorentz factor of 
the particles. Also note plus sign in front of the cyclotron term in Eq. (0). Two regimes of excitaion are 
possible: kinetic and hydrodynamic. n 

The X mode can not be excited by the Cherenkov-type interaction since it has has E _L B,v. At radii 
larger than w 50Rns (see Eq. a) the phase velocity of the of the X mode becomes smaller than the velocity 
of the primary beam with 7f, w 10^, so it can in principle be excited by cyclotron resonance, but near the 
surface of the neutron star the growth rate of the cyclotron excitation of the X mode is extremely small and 
the frequencies do not correspond to observed ones (Lyutikov 1999a). The cyclotron excitation of X mode 




^ In a hydrodynamic type instability all the particles of the beam resonate with one normal wave in the plasma, while in a 
kinetic regime only small fraction of the beam particles resonates with a given wave. 



may occur in the outer region of pulsar magnetosphere and is considered as a viable mechanism of pulsar 
radio emission generation (Kazbegi et al. 1989, Machabeli & Usov 1989, Lyutikov et al. 1999b). Thus, we 
conclude that near the surface the X mode cannot be excited and in what follows we concentrate on the 
possible excitation of the Alfven mode. 



3. Cherenkov excitation of the Alfven mode 

The Cherenkov excitation of Alfven waves has been considered in detail by Lyutikov (1999a). On a 
microphysical scale it is similar to the Cherenkov excitation of the Langmuir waves - the motion of the 
resonant particle is coupled to the component of the electric field of the wave along the magnetic field (and 
the velocity of the particle) . It is possible only for oblique propagation since Alfven waves which have a 
component of electric field parallel to the external magnetic field e^ ~ sinO. 

Using resonance condition (|^) and the low frequency asymptotics of Alfven waves (|lj) we infer that 
that the possibility of the Cherenkov excitation of the Alfven waves by the particles from the primary beam 
depends on the parameter 
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Alfven waves can be excited by Cherenkov resonance only for /i < 1, when the phase velocity of the 
fast particles may become equal to the phase velocity of obliquely propagating Alfven waves (see Figs, g 
and ||) . When /i > 1 the particles move faster that the beam and, due to the specific dispersion of Alfven 
waves, cannot resonate with them. In the outer parts of magnetosphere ( r > bORNs) parameter /i becomes 
much larger than unity /i 3> 1 so Alfven waves cannot be excited by Cherenkov interaction. 

The resonant condition for the Cherenkov excitation of Alfven waves (Eq. ^) may be solved for 
k± = ksinO: 




^^ = '^p^ 773 - rr (6) 



As expected, excitation is limited to ^ < ^^j e.i. for r < 50Rns- 

The growth rate for the Cherenkov excitation of Alfven waves (which occurs in the hydrodynamic 
regime) has been calculated in Lyutikov (1999) (see also Godfrey et al. 1975): 

_ 1 2 

_ V 3 uj^ UJGJ 3 cot 6 

where Wgj — Aire^nQj /m is the plasma frequency associated with the Goldrcich- Julian density. It is very 
similar to the growth of the Cherenkov excitation of Langmuir waves. This result is expected since the 
microphysics of the Cherenkov excitation of Langmuir and Alfven waves is the same - coupling of the 
parallel (to the magnetic field) motion of particles to the parallel component of the electric field of the 
wave. Consequently, the growth rate given by Eq. (R) for the Alfven wave excitation suffers from the same 
problem as Langmuir wave excitation: it is strongly suppressed by large Lorentz factor of the primary 
beam. We thus conclude that Cherenkov excitation of the Alfven waves is ineffective. 



4. Cyclotron excitation of the Alfven mode 

The other possibihty of excitation of Alfven waves is by the anomalous cyclotron resonance (Tsytovich 
& Kaplan 1972, Hardee & Rose 1978, Lyutikov 1999a). On a microphysical scale, during an emission at the 
anomalous cyclotron resonance, a particle undergoes a transition up in Landau levels coupling transverse 
velocity to the electric field of the wave (Ginzburg & Eidman 1959) . 

In case of low frequency {loa ^ (^p) waves, the the resonance condition for the cyclotron excitation of 
Alfven reads 

kresC cost) —^ ^ — ^ + UJEHb = 0, ^resC < UJp (8) 

If the third term is much larger than the first two (this happens for the angles of propagation larger than 
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some critical angles -^^ and F ) the resonance occurs at 

7b cos 9 sin" y I fi"^^'' ycosW sm 



kresC = I "^ '", I « ^^ ( ^—^ 1 (9) 



The condition ujres ^ ^p, which guarantees that Alfven waves are not damped by the Cherenkov interaction 
with the bulk plasma particles, requires ^ ^ cos e sin^ e ' °'' equivalently, r > 50Rns- This is a serious 
restriction: Alfven waves cannot be excited at lower altitudes since then the cyclotron resonance occurs in 
the region where Alfven waves are strongly damped due to Cherenkov resonance with the bulk particles 
(Arons & Barnard 1986). 

Cyclotron excitation of Alfven waves occurs in kinetic regime (Lyutikov 1999a). The growth rate is 

r = 5^ (10) 

where A7 is the scatter in Lorcntz factors of the resonant particles and the resonant frequency LUres follows 
from Eq. (^). Formally, growth rate for the cyclotron instability on Alfven waves (Eq (|^) is the same as 
the growth rate for the the cyclotron instability on the high frequency transverse waves (Lyutikov et al. 
1999b, Kazbegi et al. 1989), which occurs in the outer parts of the pulsar magnetosphere. The important 
difference is that the cyclotron instability on Alfven waves can occur at lower altitudes, where the density 
of the resonant particles is higher, n 

Numerically, for the emission generated at r « 50Rns and a quite narrow primary beam (A7 « 10^ - 
see Lyutikov et al. 1999a), we have 

F w lO^sec^i (11) 

Thus, we conclude that the growth rate of the cyclotron instability on the Alfven waves may be considerable 
enough to account for the high brightness pulsar radio emission. 



5. Wave conversion 

There are two fundamental problem with Alfven waves - they cannot escape from plasma and they 
are damped on the Cherenkov resonance with plasma particles when their frequency becomes comparable 



^In the case of Alfven waves we use the dispersive correction k^c^ sin^ ^/'^p instead of 'jJp/'jJ 



to LOp. As Alfven waves propagate into decreasing plasma density, their frequency will eventually become 
comparable to the plasma frequency (Barnard & Arons 1986). Before that, Alfven waves have to be 
converted into escaping modes. The conversion should take place at such radii that the local plasma 
frequency, transformed into the pulsar frame, is still larger than the observed frequencies: r < 500i?jvs- 

There are two generic types of conversion - linear and nonlinear. Linear conversion occurs in the regions 
where WKB approximation for wave propagation is not satisfied (e.g., Zheleznyakov 1996). Nonlinear 
conversion is due to the wave-wave or wave-particles interaction. In order to escape absorption at the 
Cherenkov resonance the Alfven waves should converted into either superluminous ordinary waves or into 
X modes that does not have a component of the electric field along the external magnetic field. 



5.1. Linear conversion 

Linear conversion of waves occurs when dispersion curves of two modes approach each other closely on 
the u) — k diagram. Effective conversion occurs when the distance between two dispersion curves becomes 
comparable to the "width" of the dispersion curve. Several process can contribute to the broadening of the 
dispersion curve. First, the inhomogeneity of the medium induces a width ^ ~ -^ where L is a typical 
inhomogeneity scale. Inhomogeneities can be due both to large scale (of the order of the light cylinder 
radius) density fluctuations, excited by temporal and special modulations of the flow, and due to the small 
scale plasma turbulence. Linear conversion of Alfven waves into X mode is impossible because of their 
different polarizations. Linear conversion of Alfven waves into O mode can occur only when the frequency 
of the Alfven waves approaches plasma frequency, e.i. exactly at the moment when Alfven waves become 
strongly damped at the Cherenkov resonance with the bulk plasma. The effectiveness of the conversion 
then depends on the not very well known details of the bulk plasma distribution function (e.g., distributions 
with considerable high energy tail will tend to damp the waves stronger), so no decisive argument about 
the effectiveness of such conversion can be made. 

Secondly, the wave turbulence results in effective collisions (with frequency Vc) of Alfven waves with 
each other; in this case the induced width is 5u) « Vc- □ For a given turbulent energy density Wtot the 
typical wave-wave collision frequency (based on a three wave interaction) can be estimated as 

2 



f-)^ (12) 



Effective conversion occurs when this frequency is of the order of the minimum difference between the 
dispersion curves. In the case of strong turbulence the the collision frequency becomes comparable to the 
frequency of the waves near the conversion point: v k, lo^. In that case the required wave energy density is 
approximately equal to the energy in plasma. Assuming that the energy density of plasma is of the order of 
the thermal energy, we find 

Ty*. = (^)'«nmc^ (13) 

In addition, in the pulsar plasma the energy in the secondary plasma is approximately equal to the energy 
in the beam. Then, the energy density given by Eq. ( p3| ) is of the order of the beam energy. The amount 
of energy lost by a beam due to the wave excitation depends on the complicated physics of the instability 



^Terminology here may be a bit confusing: the frequency of collisions, of course, depends on the total energy density of 
turbulence. 



saturation mechanisms, but we can reasonably expect that to be ^ 1 — 10%. In this case the the collisional 
frequency Vc ^ 0.01 — O.lujp and the collisional conversion of Alfven waves into O mode can be marginally 
effective. 

Another mechanism of linear wave transformation is related to the presence of velocity shear in the 
flow (Arons & Smith 1979, Mohajan et al. 1997, Chagelishvili et al. 1997). We can reasonably expect that 
the flow of pair plasma has some shear associated with the plasma generation in the polar gap. Given a 
strong shear Langmuir waves become coupled to the escaping O modes. Whether the shear expected in the 
polar outflow can provide enough coupling remains uncertain. 



5.2. Nonlinear conversion 

5.2.1. Wave - wave interaction 

Consider a merger of two Alfven waves into escaping X or O mode. For the three wave processes 
to take place, the participating waves should satisfy resonance conditions (conservation of energy and 
momentum) and more subtle conditions on the polarization determined by the matrix elements of the 
third order nonlinear current (Melrose 1978 Eqs. (10.105) and (10.125)). From the energy and momentum 
conservation it is easy to see that only oppositely propagating Alfven waves can merge into X or O mode. 

In calculating the matrix element of the three-wave interaction two simplifications are possible in the 
case of strongly magnetized electron-positron plasma. First, nonlinear current terms which are proportional 
to an odd power of the sign of the charge will cancel out. Secondly, we can make an expansion in powers 
I/ojb and keep the lowest order. Under these assumption the matrix element becomes 

V^-^^ (14) 

mcujB 

Given the matrix element (H) we can find the the probability of emission in the random phase approximation 
(which assumes that radiation is broadband) (Melrose 1978) 

uKV^Up = ^^^ (15) 

Then, again assuming that the energy density of plasma is of the order of the thermal energy, the 
characteristic nonlinear decay time is 

(compare with Mikhailovskii 1980). So that at the distance r « 50 — IOORns the nonlinear conversion of 
Alfven waves is only marginally possible. The nonlinear conversion of Alfven waves at lower altitudes is 
impossible because of the small growth rate. 

Thus, we come to a conclusion that though in principle the Alfven waves can be converted into 
escaping modes by nonlinear selfinteraction, the conversion must occur at comparatively large radii 
r « 50 — 100i?Ars. There is another serious problem with nonlinear conversion - it requires a presence of 
strong backward propagating waves. Presence of such backward propagating waves is not obvious in the 
pulsar magnetosphere. They should either be excited by a backward propagating fast particles, whose 




origin cannot be simply justified, or be a result of a scattering of the initial forwards propagating waves. 
The Thompson scattering is strongly suppressed in the magnetized plasma and is probably ineffective, but 
induced scattering (see below) may provide such backward propagating waves. 



5.2.2. Induced scattering 

Induced scattering of longitudinal waves in pulsar magnetosphere have been considered by Machabeli 
(1983) and Lyubarsky & Petrova (1996). They showed that it may be an effective process of transferring 
energy from the Alfven branch to the O mode. There are several problems with this mechanism. First, 
the effective induced scattering of Alfven waves on plasma particles in a superstrong magnetic field occurs 
in the same region where waves become strongly damped, the question which process is more effective 
(scattering or damping) again strongly depends on the unknown details of the plasma distribution function. 
Secondly, if scattering is effective, the induced scattering transfers energy to smallest wave vectors of the O 
mode (Langmiur condensate) . Propagation and escape of the O mode has not been properly investigated, 
but since the O mode with initially small wave vector has a very small index of refraction {n ^ 0) it will 
be strongly refracted (by the angle « tt) as it converts into a vacuum mode with (n ^ 1). This strong 
refraction would then contradict the observed narrow pulsar profiles. 

Summarizing this section, we conclude that at this point we cannot make a decisive statement whether 
linear or nonlinear conversion of Alfven wave into escaping modes is effective. 



6. Conclusion 

The results of this work confirm the conclusion of Lyutikov (1999a) that in the pulsar magnetosphere 
the electromagnetic cyclotron instabilities are the most likely candidates for the pulsar radio emission 
generation. These instabilities develop on the X and O modes (in the outer regions of the pulsar 
magnetosphere) and on the Alfven mode (possibly in the lower, r ^ 50 — IOORns, regions). Cyclotron 
instabilities on the X and O modes have smaller (than on the Alfven waves) growth rate, but generate 
waves that can directly escape from the plasma. The growth rates of the cyclotron instability on Alfven 
waves can be very large, but the complications of the wave conversion and absorption in the outflowing 
plasma, which depend on the unconstrained details of the plasma distribution function, may put the model 
based on the Alfven wave excitation at disadvantage. 

Both excitation and the nonlinear conversion of Alfven waves is possible only for r > 50 — IOORns- 
This is a comparatively large altitudes. If we suppose that the opening angle of the magnetic field lines 
controls the width of the pulsar emission beam (e.g, Rankin 1992), then the value of the opening angle 
would imply emission right near the surface r ^ Rns- If taken at a face value, this interpretation would 
exclude the Alfven waves as a source of observed radio emission. 

The comparatively large radii of emission, r > 50 — IOORns and larger, may not be that unacceptable 
from the observational point of view as well (Lyutikov 1999b). Such effects as "wide beam" geometry 
(Manchester 1995), emission bridge between some widely separated pulses, extra peaks in Crab at high 
frequencies maybe naturally explained by large altitudes of emission. The results of this work stress one 
again the difficulties of wave excitation at very small altitudes: Alfven waves, as well as ordinary and 
extraordinary modes, cannot be excited beam instabilities or converted nonlinearly into escaping modes 



at low altitudes. On the other hand, both excitation and conversion should occur at larger r > 50Rns 
altitudes which remain our prefered regions for the generation of pulsar radio emission. 

I would like to thank George Machabeli for discussions and useful comments. 
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FIGURE CAPTIONS 

Fig. 1(a). Resonances on the Alfven mode for /i < 1. 

Fig. 1(b). Resonances on the Alfven mode for /i > 1. 
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